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Abstract 

Undercooling a melt often facilitates a metastable phase to nucleate preferentially. Although the 
classical nucleation theory shows that the most critical factor for forming a metastable phase is 
the interface free energy, the crystallographic stability is also indispensable for the phase to be 
frozen at ambient temperature. In compound materials such as oxides, authors have suggested 
that the decisive factors for forming a critical nucleus are not only the free energy difference but 
also the difference of the entropy of fusion between stable and metastable phases. In the present 
study, using REFe 03 (RE: rare-earth element) as a model material, we investigate the formation 
of a metastable phase from undercooled melts with respect to the competitive nucleation and 
crystallographical stabilities of both phases. 

Introduction 


a 


A metastable phase is a phase that does not exist in thermal equilibrium state and, although 
thermodynamically unstable, can temporarily exist when some conditions are fulfilled. Research 
into the metastable phase began with Ostwald’s prediction that a phase formed first from 
supersaturated liquid is not always thermodynamically stable but is close to liquid in energy [1], 
This prediction is called “step rule”. Later on, Stranski and Totomanov [2] suggested that the 
step rule is a consequence of preferential formation of a critical nucleus of the metastable phase. 
That is, the activation energy required to form a critical nucleus, A G„ , controls the nature of the 
process. Regarding this point, the classical nucleation theory [3] states that A G n can be 
understood in terms of the interfacial free energy y between the liquid and solid phases. Turnbull 
[4] and Spaepen [5], assuming that y of a simple material such as metal is related not to the 
enthalpy change but to the entropy change at the solid-liquid interface, formulated y as 


r = 


A S f T 


(i) 


where AS/, T, Na and V,„ are the heat of fusion, the temperature of material, Avogadro number 
and the molar volume, respectively. Furthermore, Spaepen and Meyer [6] derived a, 
dimensionless solid-liquid interfacial energy, as 0.86 for fee or hep crystals and 0.71 for bee 
structures, respectively. The a-factors, which strongly depend on the structure of both solid and 
liquid phase, are to be a critical parameter to determine A G n . In fact, it has been reported that the 
phase selection of the stable y-phase or the metastable 5-phase in Fe-Ni-Cr alloys is controlled 
by a [7-9]. However, almost the metastable phases formed at the first recalescence, which 
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change into the stable phases at the second recalescence, are not frozen into ambient temperature. 
In order that the metastable phase may be frozen into ambient environment, the nucleation 
criterion of metastable phase as well as the crystallographical stability must be fulfilled. Then, in 
the present investigation, using REFe 03 as the model material, where RE means rare-earth 
elements, the phase selections not only in the nucleation stage but also the growth stage are 
discussed. 

Entropy -undercooling regime criterion of phase selection 

Before taking up the main subject of this paper, we mention again the hypothesis that, in ionic 
crystals, AS/ is to be a dominant factor in the determination of y. 

Spaepen [10] and Granasy [11], almost at the same time, developed rather similar models that 
y at equilibrium state is given by 

y= J v (H(r)-TS(r)y r , (2) 

where H(r) and S(r) are cross-interfacial enthalpy and entropy. Figure 1 shows a schematical 
illustration of Eq. 2, where (a) shows the change of the atomic order and of the order parameter 
in solid and liquid, and (b) schematic representation of the change in the enthalpy H(r) and the 
□ product T E S(r ) of melting temperature and entropy at the solid-liquid interface. The colored area 
corresponds to the approximate value of the interfacial energy, showing that a is not a 
dimensionless interfacial energy but a dimensionless interface thickness. This model can 
qualitatively infer the temperature dependence of the solid-liquid interfacial free energy. 
However, we need analytical or numerical forms of H(r ) and S(r) in order to evaluate the 
interface thickness that is the key parameter of the model. 

On the other hand, in a material having faceted interface, the order parameter representing the 
regularity of the atomic arrangement rapidly changes as the interface is crossed as shown in Fig. 
2. In this case, H{r) can be approxiomated by the near- step function as shown in Fig. 2(b). 
Consequently, the interfacial energy is approximated with a triangle if we assume S(r ) as a linear 
function of r at the cross-interface region. This result also means that a is a dimensionless 
interface thickness rather than a dimensionless interfacial energy. The recent numerical 
calculation of a based on the molecular dynamics and the density functional analysis suggests 



Figure 1 . Schematic representation of the nonfaceted interface between solid and liquid, (a) 
change of the atomic order and of the order parameter in solid and liquid, (b) schematic 
representation of the change in the enthalpy H{r) and the product T E S(r ) of melting 
temperature and entropy at the solid-liquid interface. The colored area corresponds to the 
approximate value of the interfacial energy, showing that a is not a dimensionless interfacial 
energy but a dimensionless interface thickness. 
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Figure 2. Schematic representation of the faceted interface between solid and liquid, (a) 
change of the atomic order and of the order parameter in solid and liquid, (b) schematic 
representation of the change in H(r ) and T E S(r). The colored area corresponds to the 
approximate value of the interfacial energy, suggesting the interface thickness is 
approximately half of the atomic layer spacing. 


Figure 3. Schematical image of temperature dependency of free energy in liquid and solid 
phases (stable and metastable phases). The entropies of three phase are related Si>S,„,>S s . 
Therefore, as for entropy of fusion AS/ caused by solidification, the relation AS/P>AS/ ms 
becomes valid. 

the interface thickness is approximately half of the atomic layer spacing [12], which implies that 
the model shown in Fig. 2 is qualitatively valid even in the nonfaceted interface. 

Figure 3 schematically shows the thermodynamic relation when considering the step rule [13]. 
The figure depicts temperature and free energy of both liquid and solid phases (stable and 
metastable phases). The reason why liquid phase changes to solid phase is that the free energy of 
the liquid phase becomes larger than that of solid phase. The energy-balance point of both phases 
is the melting point. Comparing the metastable phase to the stable solid phase in terms of free 
energy, the free energy of the metastable phase, G ms , is larger than that of stable phase, G s , (the 
subscripts s and ms mean stable phase and metastable phase, respectively). Therefore, the melting 
point of the metastable phase, T E , ms , becomes lower than that of stable phase, T E , S . Meanwhile, 



T T 

Temperature Ems Es 
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the absolute value of the gradient of each curve (temperature coefficient of Gibbs free energy) in 
Fig. 2 corresponds to the entropy when pressure is constant. From the figure, we can see the 
relation S/>S mi >S s between entropies of liquid phase Sl, stable phase S s and metastable phase S ms . 
Therefore, for the change of entropy, AS/, caused by melting, we can find a relation A Sf s >A Sf ms 
(i.e., the change of the entropy is smaller when the liquid phase changes to metastable phase). 
From the relative relations of the three phases above, we can see that the metastable phase is to 
be a higher entropy phase than the stable phase [14], 

Factors determining the entropy of material are first, density of material and secondly, 
symmetry of arrangement of atoms and/or molecules making up the material. Therefore, high 
entropy phase is liquid rather than solid, and gas rather than liquid. Among solid phases, it is 
guessed that the low-density phase becomes higher-entropy phase. In conclusion, we can say that 
the metastable phase is lower density, higher symmetric material than the stable state. 


Experimental procedure and results 

Spherical samples of REFe0 3 were prepared from high purity (99.99%) RE2O3 and Fe203 
powders. Levitation and melting of samples were carried out by an aerodynamic levitator, ADL, 
which was designed in order to solidify undercooled melts under the precisely controlled P 02 . 
Details of the sample preparation and experimental facility are shown elsewhere [15, 16]. 

Goldschmidt [17] discussed the stability of the perovskite (ABO3) structure using the 
tolerance factor, TF: 

R + R 

TF =-j=r± 2 — . (3) 

V 2(R b + R 0 

In the present investigation, ionic radii of Ra, Rb and Rq correspond to those of rare-earth 
element, iron and oxygen, respectively. From the systematic investigation, he summarized that 
the perovskite structure is stable at TF > 0.8 and contrary unstable at TF < 0.8. Using Shannon 
ionic radii from La (0.1216 nm) to Lu (0.1032 nm), TF ’ s for the REFe03 system were calculated 
to be 0.905 for LaFe03 to 0.841 for LuFeCU Therefore, the perovskite structure is expected to be 
stable in the REFe03 system. 

Figure 4 shows SEM micrographs of samples of REFe0 3 solidified in containerless conditions. 
As shown in the photographs, their surface profiles vary according to the type of rare-earth 
elements. The surface of LaFe03 is nonfacetted and spherical while that of LuFe0 3 is facetted and 
polyhedral [18]. Note that the different surface features result from differences in crystal structure, 
not differences in rare-earth element. Specifically, the lattice structure of the LaFe03 sample is 
orthorhombic, the space group of which is Pbnm (o-REFe0 3 ), while the LuFe03 sample is a 
hexagonal- symmetric of P6 3 cm (/t-REFe0 3 ). As the /t-REFe0 3 phase has a 10%~20% smaller 



Figure 4. Surface profiles of REFe0 3 (RE=La, Lu) formed by containerless process. LaFe0 3 
with large ionic radius has a smooth and spherical surface while LuFe0 3 with small radius 
has a rugged and polyhedral surface. 
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density than that of perovskite [19], it is estimated that the hexagonal crystal is a higher entropy 
phase than perovskite. In other words, the hexagonal crystal should have intrinsically become 
stable perovskite. However, having been largely undercooled to below Te , ms indicated in Fig. 2 by 
the containerless process, hexagonal crystals of high-entropy phase grew as metastable phase 
(Fig. 5). In fact, when we forced it to solidify at a temperature of around Te, s even by the same 
containerless process, stable-phase perovskite appears. 

Figure 5 shows typical images taken sequentially during recalescence in samples of LaFeC> 3 , 
GdFe0 3 and YFe0 3 , each of which is processed at oxygen environment. The elapsed time 
indicated in each image was set to 0 s for one frame before the nucleation. At oxygen 
environment, although single recalescence that can be ascribed to the phase transition from 
undercooled melt to equilibrium perovskite phase was observed in LaFeC >3 samples, double 
recalescences were observed in GdFe0 3 and YFe 03 samples, where a primary phase was 
solidified from the undercooled melt and then the secondary phase with higher brightness was 
initiated at the interface between melt and the primary phase. The high brightness implies that 
the melting temperature of the secondary phase was much higher than that of the primary phase. 
This result indicates that decrease of TF facilitates the undercooled melt to solidify into the 
metastable /?-REFe0 3 phase rather than the stable o-REFeCF phase. According to this result, 
reduction of the oxygen partial pressure P 02 is expected to extend the range of TF for metastable 
h- REFe 03 phase to be formed, because decreasing P 02 increases the amount of Fe 2+ (0.078 nm 
for CN=6) with larger ionic radius than that of Fe 3+ (0.0645 nm for CN=6). 


0 s 0.001 0.002 0.005 



0 s 0.019 0.027 0.038 0.049 



0 s 0.020 0.026 0.051 0.083 



Figure 5. Sequence photographs of HSV images taken during recalescences in the REFe0 3 (R = La, 
Gd and Y) samples processed at oxygen environment. Although in LaFe0 3 single recalescence was 
observed, in GdFe0 3 and YFe0 3 double recalescences indicating the formation of metastable phases 
were observed. 

In addition, YbFe 03 was used to study the effect of ionic radii of RE elements on the 
formation of metastable phases, because the ionic radius of Yb 3+ (0.1042 nm) was slightly larger 
than that of Lu 3+ (0.1032 nm). Figure 6 shows the XRD patterns of YbFe 03 samples processed at 
controlled P 02 . At 10 5 Pa of P 02 , the stable orthorhombic phase (o-YbFe0 3 ) was formed at the 
second recalescence as in the cases of GdFe 03 and YFeC> 3 . At 10 4 Pa of P 02 , however, the 
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metastable hexagonal phase (/j-REFeCF) remained, forming the dual phase with o-YbFeCF, and 
at 9xl0 3 Pa, the o-YbFeCF phase did not appear [20]. 

These results suggest that the decrease of Po 2 facilitates the undercooled melt to solidify to 
metastable /r-REFeCF phase rather than the stable o-REFeCF phase, particularly in samples with 
RE 3+ of relatively small ionic radius. 

Figure 7 shows the relation between the recalescence results and P 02 as a function of the ionic 
radii of RE 3+ , in which Shannon ionic radii for CN=9 were used. Decrease of Poj extends the 



20 / degrees 

Figure 6. XRD patterns of the YbFe0 3 samples processed at Po 2 =10 5 Pa, Po 2 =10 4 Pa and 
Po 2 =9x10 3 Pa, respectively. At 10 5 Pa of Po 2 , the stable orthorhombic phase (o-YbFe0 3 ) was 
formed at the second recalescence as in the cases of GdFe0 3 and YFe0 3 . At 10 4 Pa of Po 2 , 
however, the metastable hexagonal phase (/z-REFe0 3 ) remaind, forming the dual phase with o- 
YbFe0 3 , and at 9xl0 3 Pa, the <?-YbFe0 3 phase did not appear. 



Atomic Number 
57 


62 63 64 


68 


70 71 


1.22 


1.14 


1.02 


1.18 


1.10 


1.06 


Ionic Radius of R 3+ (A) 


Figure 7. Relation between the recalescence results and Po 2 as a function of the ionic radii 
of RE 3+ , in which Shannon ionic radii for CN=9 were used. 
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range of TF for metastable /i-REFeCF phase to be formed. 


Discussion 


Bertaut et al. [21] and Yakel et al. [22] have first reported the two hexagonal modifications in 
the ABO3 systems, the space groups of which are PG^hnmc and PG^cm, respectively. In the h- 
REMn03 system, the PG^cm type modification was formed as a low temperature phase for RE 3+ 
with small ionic radius (Ho-Lu, Y and Sc), whereas the PG^/mmc type modification was reported 
as a high temperature phase. In our experiment, although the space group of the metastable h- 
REFe03 phase belonged to PG^cm, the high temperature phase can be deduced to belong to the 
PG^mmc space group because the ionic radius of Fe 3+ is as same as that of Mn 3+ (0.0645 nm for 
6 coordination) [23]. Hence, in this investigation, the geometrical analysis of the atomic 
configuration in /i-RFFe0 3 is developed on the assumption that the space group of the primary 
phase is PG^/mmc. 

The atomic configuration of the PG^cm modification in ABO3 system can be described as a 
dense oxygen-ion packing (ABCACB) with B 3+ ions having coordination number CN=5 (five- 
fold distorted trigonal bipyramidal coordination), and A 3+ with CN=7 (seven-fold monocapped 
octahedral coordination), forming a noncentrosymmetric structure. On the other hand, 
centrosymmetric PG^/mmc is assumed to be described simply with B 3+ ions of undistorted CN=5 
and A 3+ of CN=6 (octahedral coordination). 

Figure 8 shows the geometrical configuration among RE 3+ (CN=6), Fe 3+ (CN=5), and O 2 ' in a 
space group of PG^/mmc, in which the constituent ions are packed without any spacing between 
neighboring ions. As shown in this figure, the relation among the ionic radii of constituent ions 
of /t-REFe0 3 is expressed as 

4g , 

R +R = (R + R\ (4) 

Therefore, the /r-REFeCF phase will be ideal when the next equation is fulfilled, 

0.87 (5) 

results showed that the /?-GdFe0 3 phase is formed because TF of which is 
h- EuFe03 is not formed because of slightly high TF of 0.871. This suggests 


TF = 


F 


H 


The experimental 
0.867. Whereas, the 


□ 



Ke + *0 = ~(^Fe + Ao) 

^±Aq V3 
V2 (R b + Ro) 2 


Figure 8. Geometrical configuration among RE 3+ (CN=6), Fe 3+ (CN=5), and O 2 ' in a space 
group of I’hTirwic. in which the constituent ions are packed with no space between 
neighboring ions. The /i-REFe0 3 phase will be the ideal at TF = 0.87. 
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that the aforementioned condition expressed by Eq. 5 is to be the criterion for the metastable h- 
REFeOs phase to be formed. 


Conclusion 

Using REFe 03 (RE: rare-earth element) as a model material, containerless solidification for 
forming a metastable phase from undercooled melts was carried out as a function of P 02 (oxygen 
partial pressure). Based on the geometrical consideration on the ionic radii of constituent ion, 
RE 3+ , Fe 3+ and O 2 ", tit was derived that the criterion for metastable hexagonal phase is also 
expressed by tolerance factor, TF: 

TF < 0.87. 

Experimental result well agreed with this criterion under reduced P 02 as well as for ambient 
conditions. 
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